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ABSTRACT 

We demonstrate that a log-linear relation does not provide an adequate description 
of the correlation between the masses of Super-Massive Black-Holes (SMBH, Mbh) 
and the velocity dispersions of their host spheroid (cr). An unknown relation between 
log Mbh and log a may be expanded to second order to obtain a log-quadratic relation 
of the form log (Mbh) = a + /31og (tr/200kms _1 )-(-/32[log (er/200kms -1 )] 2 . We perform 
a Bayesian analysis using the Local sample described in Tremaine et al. (2002), and 
solve for /3, /?2 and a, in addition to the intrinsic scatter (6). We find unbiased param¬ 
eter estimates of f3 = 4.2 ± 0.37, /?2 = 1-6 ± 1.3 and 6 = 0.275 ± 0.05. At the 80% level 
the Mbh — cr relation does not follow a uniform power-law. Indeed, over the velocity 
range 70km/s^ a ^380km/s the logarithmic slope d log Mbh /d log a of the best fit 
relation varies between 2.7 and 5.1, which should be compared with a power-law esti¬ 
mate of 4.02 ± 0.33. The addition of the 14 galaxies with reverberation SMBH masses 
and measured velocity dispersions (Onken et al. 2004) to the Local SMBH sample 
leads to a log-quadratic relation with the same best fit as the Local sample. Further¬ 
more, assuming no systematic offset, single epoch virial SMBH masses estimated for 
AGN (Barth et al. 2005) follow the same log-quadratic Mbh — cr relation as the Local 
sample, but extend it downward in mass by an order of magnitude. The log-quadratic 
term in the Mbh — cr relation has a significant effect on estimates of the local SMBH 
mass function at Mbh 1O 9 M 0 , leading to densities of SMBHs with Mbh ^ 1O 1O M 0 
that are several orders of magnitude larger than inferred from a log-linear Mbh — cr 
relation. We also estimate unbiased parameters for the SMBH-bulge mass relation 
using the sample assembled by Haering & Rix (2004). With a parameterization 
log (Mbh) = Ckbuige + /3buige log ( M hu i ge /l0 n M & ) + /3 2 ,bulge [log (Mb u ig e /1O 1:L M 0 )] 2 , we 
find /3bui g e = 1-15 ± 0.18 and /? 2 ,bulge = 0.12 ± 0.14. We determined an intrinsic scatter 
<5buige = 0.41 ± 0.07 which is ~50% larger than the scatter in the Mbh — cr relation. 
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1 INTRODUCTION 

Observations of bulges in nearby galaxies reveal the pres¬ 
ence of massive dark objects whose dynamical influence on 
the surrounding stars is consistent with their being SMBHs 
(e.g. Kormandy & Richstone 1995). Moreover, the masses of 
these SMBHs correlate with properties of the host galaxy, 
including the luminosity of the bulge (Kormendy & Rich- 
stone 1995), the mass of the bulge (Magorrian et al. 1998), 
the stellar velocity dispersion of the bulge (Ferrarese & Mer¬ 
ritt 2000; Gebhardt et al. 2000) and the concentration of the 
bulge (Graham et al. 2002). The tightest relation, with in¬ 
trinsic scatter of ~ 0.3 dex, appears to be between SMBH 
mass and bulge velocity dispersion (Tremaine et al. 2002, 
hereafter T02). Significant effort has been invested in de¬ 


termining parameters that describe this relation, which is 
usually parameterised using the log-linear functional form 

log(Mbh) = a + /Hog(<j/200kins _1 ). (1) 

The value of the power-law slope /3 has been a matter of 
some debate (see T02 for a review). Recent estimates have 
j3 = 4.02 ± 0.33 (T02) and j3 = 4.83 ± 0.43 (Ferrarese & 
Ford 2004). These values differ by ~ 2-sigma, a difference 
which may be attributable to systematic differences in the 
velocity dispersions used by different groups (T02). 

Although the connection between SMBHs and their 
host galaxies is not yet clear, it seems very likely that their 
evolution is closely intertwined. It also seems likely, given the 
small intrinsic scatter in the relation that the value of /3 will 
yield important clues regarding the physics of SMBH evolu- 
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tion (e.g. Silk & Rees 1998; Wyithe & Loeb 2003; King 2004; 
Miralda-Escude 2004; Saznov et al. 2005). Moreover, studies 
of SMBH demographics (e.g. Yu & Tremaine 2002; Shankar 
et al. 2004) rely on this relation to estimate quantities of 
astrophysical interest like the local SMBH mass density. For 
these reasons, and given the heroic efforts that have been 
made to measure SMBH masses and host velocity disper¬ 
sions, it is important to make statistically robust estimates 
of parameters that describe the Mbh — cr relation. 

Rather than assume a powerlaw relation between Mbh 
and a we take a more general approach in this paper. Sup¬ 
pose we have an unknown relation log(Afbh) = /(log(<r)). 
This relation may be expanded in a Taylor series in log a 
about a = 200km/s, which yields to second order 

log(Afbh) = a + /Jlog(a-/200kms -1 ) 

+ fo [log(a/200kms -1 )] 2 . (2) 

The zeroth and first order coefficients in this expansion cor¬ 
respond to a and (3 in the usual log-linear relation. However 
in addition to a and f3 we also consider the possibility of 
a log-quadratic term. The parameter /?2 therefore provides 
a measure of the deviation from a pure powerlaw of the 
Mbh — cr relation. 

We show that there is a non-linear contribution to the 
log Mbh — log a relation, and that the assumption of a log- 
linear relation has biased estimates of the power-law slope. 
For the Local sample of T02 we show that the variation in 
log-linear slope over the range of a in the sample is larger 
than both the statistical uncertainty in the slope of the log- 
linear relation and the 2-sigma difference between the es¬ 
timates of T02 and Ferrarese (2002). Note that (with the 
exceptions of § 5 and the appendices) we restrict our atten¬ 
tion to the statistical bias in estimation of the parameters 
(3 and @2 from SMBH masses and galaxy velocity disper¬ 
sions summarised in T02. For further discussion of possible 
systematic errors in the measured parameters Mbh and cr 
themselves we refer the reader to T02 and to Merritt & Fer¬ 
rarese (2001). 

An outline of the paper is as follows. We first discuss 
possible bias in the local sample of SMBHs with respect to 
parameter fitting (§ 2). Then in § 3 we repeat the parameter 
fitting analysis of T02 and demonstrate that the value of f3 
in a log-linear relation depends systematically on which low 
velocity dispersion galaxies are included in the sample. We 
then show that a log-quadratic form for the Afbh — er relation 
provides an improved fit which is not systematically sensitive 
to the galaxy sample (§ 4). The dependence of this conclu¬ 
sion on the method of SMBH mass estimation and resolution 
of the SMBH sphere of influence, as well as the definition of 
the velocity dispersion is discussed in § 5. A Bayesian ap¬ 
proach to parameter estimation in the Mbh — cr relation and 
quantification of the systematic bias in parameter estima¬ 
tion are described in § 6. This approach allows an unbiased 
analysis for the Local sample which is presented in § 7. In 
§ 8 we apply our analysis to the relation between SMBH 
and bulge mass using the sample of Haering & Rix (2004). 
We present some discussion in § 9, before summarising our 
conclusions in § 10. Finally, in two appendices we analyse 
alternative samples to the one described in T02. 


2 PARAMETER ESTIMATION IN THE 

LOCAL SMBH SAMPLE 

A first question that should be considered concerns whether 
or not the local SMBH sample, and therefore the result¬ 
ing estimates of the Mbh — cr relation are statistically fair 
with respect to parameter estimation. The best-fit log- 
linear Mbh — cr relation (T02; Ferrarese 2002) drops below 
~ 2 x 1O 6 M0 for velocity dispersions smaller than ~ 70km/s, 
still within the range of a in the Local sample. However no 
kinematic detections of SMBH masses have been published 
below 2x 10 6 Mq. Does this lack of low mass SMBHs result in 
biased estimates of the slope in a log-linear Mbh — cr relation? 
Ferrarese (2003) has suggested that the reason for the dearth 
of SMBHs with M < 10 6 Mq is partly instrumental sensitiv¬ 
ity. SMBHs with masses smaller than 10 6 Mq are expected 
to reside in dwarf galaxies, with velocity dispersions below 
50km/s. The sphere of influence of these SMBHs would be 
observable only in the most nearby cases. Ferrarese (2003) 
has produced a scatter plot of the local SMBH population 
in SMBH mass and distance based on a combination of the 
CfA redshift survey (Huchra et al. 1990) and the relation 
between SMBH mass and bulge luminosity. For the nearest 
group of galaxies at ~ 3Mpc, resolving the SMBH sphere 
of influence limits detection to SMBHs with masses in ex¬ 
cess of 2 x 10 6 Mq. Ferrarese (2003) goes on to suggest that 
in the nearby dwarf galaxies, individual stars become re¬ 
solved, but are too faint to allow dynamical studies, and 
that enlargement of the kinematically detected sample to in¬ 
clude SMBHs below ~ 1O 5 A70 will require campaigns using 
telescopes with apertures in excess of 8m, combined with 
resolution of better than 0.02”. One might therefore sup¬ 
pose that there is a selection bias in the SMBH sample for 
SMBHs that are larger than ~ 10 6 Mq. Such a selection bias 
would lead to biased estimates of the Mbh — cr relation. How¬ 
ever inspection of the Local sample shows that one of the 
smallest SMBHs observed (N7457) lies at a distance well 
beyond where its sphere of influence would have been re¬ 
solved by HST (Ferrarese 2003). Moreover, in the velocity 
range 70km s -1 < <r <380kms _1 there are no galaxies in 
which a kinematic search has been conducted by the Nuker 
team, and where a SMBH was not found (S. Tremaine pri¬ 
vate communication). Since these galaxies are not selected 
based on their expected SMBH mass, there should be no 
bias against selection of low mass (or high mass) SMBHs at 
a fixed velocity dispersion. One can therefore estimate the 
expected SMBH mass given a velocity dispersion from the 
local SMBH sample, without bias introduced through selec¬ 
tion. 

On the other hand, the sample is biased if one wants 
to estimate the velocity dispersion at a fixed SMBH mass. 
The low surface brightness of galaxies with large velocity 
dispersions results in a maximum velocity dispersion where 
SMBHs can be measured (Ferrarese 2003). This combined 
with the exponential decline in the number density of such 
galaxies may result in an upper cutoff in the velocity disper¬ 
sions found in the observed sample. A parametric fit that 
treated velocity dispersion rather than SMBH mass as the 
dependent variable would therefore lead to a biased estimate 
of the parameters in the Mbh — cr relation. 

Clearly the selection function for the sample of SMBHs 
is complex. In particular, the distribution of velocity disper- 
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Figure 1. Left: The Albh — cr relation from Tremaine et al. (2002). Right: The 16th, 50th (solid line) and 84th percentiles of the probability 
distribution for 0 obtained by minimising y 2 as a function of the minimum a considered in the sample. The best fit relations for these 
minimum dispersions are shown in the left panel, and demonstrate the effect of the lowest velocity points on the most likely fit. 


sions in the sample does not reflect the overall distribution 
of galaxies. The non-uniformity of the distribution of veloc¬ 
ity dispersions results in statistical inferences for parameters 
in equation (1) that vary with the statistic used (Merritt & 
Ferrarese 2001). This implies that some or all of the esti¬ 
mates are biased. In order to make an unbiased estimate 
for a parameter like the power-law slope 0 one must de¬ 
termine the bias inherent in the particular statistic used. 
In this paper we overcome the selection bias introduced by 
the inhomogeneous distribution of velocity dispersions via 
an investigation of the SMBH mass vs. velocity dispersion 
relation in mock Monte-Carlo samples. 

Given that we have access to a sample which is unbiased 
in Mbh as a function of a (the Local sample, summarised in 
T02), and in the absence of a theory where velocity disper¬ 
sion is regulated by SMBH mass (theoretical prejudice), we 
follow tradition and restrict our attention to parameter fit¬ 
ting where log Mbh is the dependent variable as a function 
of log a. 


3 THE LOG-LINEAR M BH - cr RELATION 

Tremaine et al. (2002) have compiled a list of N g = 31 galax¬ 
ies (the Local sample) with reliable determinations of both 
SMBH mass (Mbh) and central velocity dispersion (cr, de¬ 
fined as the luminosity weighted dispersion in a slit aperture 
of half length R e , the effective radius of the spheroid). The 


sample* is shown in the left panel of figure 1 which illus¬ 
trates the correlation between Mbh and cr. 

We begin by repeating the analysis of T02 who estimate 
the parameters a and 0 through minimisation of a \ 2 vari¬ 
able that accounts for uncertainties in both a and Mbh- The 
\ 2 variable used is 

N g 

2 _ v-' {yi-a~ 0Xj) . 

h e2 vi+^ i ’ 

where yr and Xi are the logarithm of SMBH mass in so¬ 
lar masses and the logarithm of velocity dispersion in units 
of 200km/s respectively. The variables e X i and e y i are the 
uncertainties in dex for these parameters. This expression 
(equation 3) is symmetric in log Mbh and log cr. One might 
expect this to be a favourable property since the fit does not 
include any preconceived notions of the physical origin of the 
relation. In T02 an estimate of the intrinsic scatter (tintrins) 
in the M b h — cr relation [as defined in the y (or ALbh) direc¬ 
tion] was established by adding ef ntrins to the denominator in 
equation (3). The intrinsic scatter that resulted in a reduced 
minimum \' 2 of unity was tintrins ~ 0.27 dex. The best fit 
solution for a linear Mbh — cr relation from T02 has a = 8.13 
and 0 = 4.02, resulting in residuals for the three smallest 
galaxies that are greater than zero. In addition, the largest 
three galaxies also have residuals that are greater than zero. 


* In this paper we refer to the Local sample as the sample 
of SMBHs and effective velocity dispersion (with uncertainties) 
listed in table 1 of T02. The exception is the Milky-Way galaxy 
for which we use the updated estimate for SMBH mass of Mbh = 
(3.7 ± 1.5) X 10® Mg (Schodel et al. 2002). References for SMBH 
masses in table 1 of T02 are Verolme et al. (2002); Tremaine 
(1995); Kormendy & Bender (1999); Bacon et al. (2001); Geb- 
hardt et al. (2003); Bower et al. (2001); Greenhill & Gwinn (1997); 
Sarzi et al. (2001); Kormendy et al. (1996); Barth et al. (2001); 
Kormendy et al. (1998); Gebhardt et al. (2000); Herrnstein et al. 
(1999); Ferrarese, Ford, & Jaffe (1996); Cretton & van den Bosch 
(1999); Harms et al. (1994); Macchetto et al. (1997); Ferrarese & 
Ford (1999); van der Marel &; van den Bosch (1998); Cappellari 
et al. (2002). 
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This behaviour is symptomatic of a scenario where a linear 
relation has been fitted to a non-linear sample. 

The sample described in T02 is dominated by galax¬ 
ies with velocity dispersions in the range ~ 130 — 250km/s. 
There are a handful of SMBHs at the centers of galaxies 
with smaller velocity dispersions, including the Milky-Way. 
These SMBHs all have masses in excess of Mbh ~ 2x 10 6 Mq. 
The four galaxies with the lowest dispersions have a large 
influence over the slope inferred for a log-linear relation. Fig¬ 
ure 1 shows how the estimate of (3 varies as the galaxies with 
the lowest velocity dispersions are removed from the sample. 
Two trends are apparent. Firstly, the uncertainty in (3 is re¬ 
duced by more than 50% by the presence of the smallest few 
galaxies. Secondly, the lowest velocity dispersions reduce the 
estimate of the slope from (3 ~ 4.5 to (3 ~ 4. The estimate 
remains near (3 ~ 4.5 ± 0.5 as the next 6 smallest veloci¬ 
ties are removed from the sample. The systematic variation 
may be seen visually in the left panel of figure 1 , where the 
best-fit relations are plotted over the corresponding ranges 
of velocity dispersion. 


4 A LOG-QUADRATIC M BH a RELATION 

The systematic trend of a slope that increases by 1.5-sigma 
when the smallest 3 of N g = 31 galaxies are removed from 
the sample points to a fit that is systematically biased. One 
possibility is that a log-linear relation is not a good descrip¬ 
tion of the data set. Note that a poor fit cannot be identified 
through a reduced value of x 2 that is significantly greater 
than unity. This is because the value of intrinsic scatter 
(tintrins) is determined by the condition that the reduced 
value of x 2 equal unity for the best fit solution. However 
we can try to find a functional form that provides a bet¬ 
ter description of the data. One possible scenario involves a 
characteristic velocity dispersion, with different powerlaws 
for large and small galaxies. To test this idea we attempt to 
locate a value of a around which the slope of the relation 
changes by fitting two power-law slopes, one either side of a 
characteristic velocity dispersion. We find that the fit is im¬ 
proved by allowing the slope at high a to be steeper than at 
low a, and that the data prefers a break velocity somewhere 
in the range 150 —350km/s (1-sigma). Of course, while there 
is at least some theoretical motivation for a log-linear rela¬ 
tion, there is no reason to suppose that the A/bh — a relation 
should be a double power-law. 

In the following we take a more general approach. As 
discussed in the introduction, a general relation log Mbh = 
/(logo-) can be expanded in a Taylor series to second order 
yielding 

log(Mbh) = a + /31og(o-/200kms _1 ) 

+ (32 [log(o-/200kms _1 )]“ . (4) 

Here the coefficients (3 and (32 represent log-linear and log- 
quadratic contributions to the Mbh — cr relation. We have 
(arbitrarily) expanded about a = 200 km/s, which corre¬ 
sponds to the median velocity dispersion in the relation. The 
value of p2 provides a measure of whether or not a log-linear 
relation provides a good description of the data. 

We have repeated the x ' 2 minimisation using equa¬ 


tion (4). The 1, 1.5, 2 and 2.5-sigma ellipsoids^ for (3 and (32 
are shown in the upper right-hand panel of figure 2 for the 
Local sample. An intrinsic scatter of emtrins = 0.25 dex re¬ 
sults in a reduced x 2 of unity for the best fit solution. Thus 
a log-quadratic fit admits a slightly smaller intrinsic scat¬ 
ter than the log-linear relation. The solution shows evidence 
for a positive log-quadratic term at around the 85% level. 
The most likely solution has (3 = 4.3 and (32 = 1-9 with 
1-sigma uncertainties of 0.35 and 1.1 respectively. The best 
fit log-quadratic relation is plotted in the upper-left panel 
of figure 2 . 

Parameter estimation using a log-linear relation corre¬ 
sponds to the conditional probability for (3 given (32 = 0. 
The horizontal grey line in the upper right panel of figure 2 
shows the cut through the bi-variate probability distribution 
for (3 and (32- The contours of the bi-variate distribution 
cross this line centered around (3 = 4, which is consistent 
with expectations from the log-linear fit. However the most 
probable value (3 = 4.02 for the log-linear fit lies near the 1.5- 
sigma contour of the log-quadratic fit. This illustrates the 
point that a log-quadratic form provides a much improved 
description of the Local sample. 

Unlike the log-linear case, the four galaxies with the 
lowest dispersions do not have a large influence on the log- 
quadratic fit. Figure 2 shows that the estimates of (3 and (32 
do not vary systematically as the galaxies with the lowest 
velocity dispersions are removed from the sample. The most 
likely solution is similar for samples including all galaxies 
and for samples including only the largest 20 galaxies, in¬ 
dicating that there is evidence for a log-quadratic term in 
the main group of massive galaxies, and that the deviation 
from a power-law is not dominated by inclusion of low mass 
galaxies in the sample. The variation of the best-fit relation 
as small galaxies are removed from the sample may be seen 
visually in the left panel of figure 2 , where the best-fit rela¬ 
tions are plotted over the corresponding ranges of velocity 
dispersion. 

The above discussion relates to an expansion of the gen¬ 
eral relation log Mbh = /(log a) which is truncated at second 
order. We found that the data prefers a non-zero contribu¬ 
tion from a quadratic term. Before continuing we mention 
the possibility of a non-zero contribution from an additional 
log-cubic term. We have repeated the above analysis for the 
third order relation 

log(Afbh) = a + /31og(a-/200kms _1 ) 

+ 02 [log(a/ 200 kms -1 )] 2 

+ (3z [log(a/200kms -1 )] 3 . (5) 

By minimising x 2 over the other parameters, we find a so¬ 
lution /3 = 4.33 ± 0.6, (32 = 1.85 ± 2.0 and jh = 0.3 ± 6.0. 
Not surprisingly there is a large correlation between the odd 
terms (3 and /J 3 . This correlation increases the uncertainties 
on individual parameters. Despite having an additional free 

t These ellipsoids correspond to loci where the difference between 
X 2 (/3, 02) and the minimum x 2 is 1, 2.71. 4 and 6.63 respectively. 
When these contours are projected onto an individual parameter 
axis, they correspond to the 68%, 90%, 95% and 99% confidence 
intervals on those parameters. In the text we refer to these as the 
1. 1.5, 2 and 2.5-sigma error elipsoids. 
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Figure 2. x 2 minimisation of the log-quadratic M\ ) } 1 — a relation. Upper Right: Contours of x 2 (minimised over a ) of 0 and 0 2 - When 
projected, the contours correspond to the 1, 1.5, 2 and 2.5-sigma uncertainties on individual parameters. The point shows the most likely 
solution. Lower Left: The 16th, 50th (solid line) and 84th percentiles of the probability distribution for 0 obtained by minimising x 2 as 
a function of the minimum a considered in the sample. Lower Right: The 16th, 50th (solid line) and 84th percentiles of the probability 
distribution for 02 obtained by minimising x 2 as a function of the minimum a considered in the sample. The best fit relations for these 
minimum dispersions are shown in the upper-left panel, and demonstrate the effect of the lowest velocity points on the most likely fit. 


parameter, the log-cubic relation requires a larger intrinsic 
scatter to achieve a \' 2 of unity than the log-quadratic rela¬ 
tion. Moreover there is no evidence from the x 2 minimiza¬ 
tion for a log-cubic term in the Local sample. Therefore in 
the remainder of the paper we restrict our attention to the 
log-quadratic form for the Mbh — cr relation. 


5 CHOICE OF SMBH SAMPLE AND 
DEFINITION OF a 

The Local SMBH sample described in T02 contains galaxies 
with SMBH mass estimates obtained from data that did not 
resolve the sphere of influence (e.g. Marconi & Hunt 2003), 
and the accuracy of these masses has been called into ques¬ 
tion (e.g. Merritt & Ferrarese 2001; Ferrarese & Ford 2004; 
Marconi & Hunt 2003). In addition, the sample has SMBHs 
with masses inferred using four separate techniques, includ¬ 
ing stellar dynamics, stellar proper motions, astrophysical 
masers and dynamics of gaseous disks. Moreover, adding to 


this inhomogeneity are the different definitions for the veloc¬ 
ity dispersion that have been supported by different groups. 
In this section we investigate the effect that these properties 
of the sample may have on the conclusion of a log-quadratic 
Mhh — cr relation. These issues are also discussed further in 
two appendices (§ A and B). 

Some of the original SMBH mass estimates based on 
ground based stellar dynamical data (Magorrian et al. 1998) 
lie systematically above the Mbh — cr relation defined by 
other galaxies in the local sample with more secure mass 
estimates by as much as 2 orders of magnitude (Ferrarese 
& Merritt 2000). The correlation between the offset and the 
galaxies distance has been interpreted as evidence that the 
lack of resolution in these early observations resulted in a 
systematic resolution dependent error in the Magorrian et 
al. (1998) masses (Merritt & Ferrarese 2001). Of course it 
is not clear that one can assume masses to be incorrect be¬ 
cause they don’t agree with the log-linear Mbh — cr relation 
as described by the remaining SMBHs. This would require 
the argument that the observational mass estimates must 


© 2005 RAS, MNRAS 000, 1-21 













6 


Wyithe 



a (km/s) 



Other 


Stellar 



(S 


(N 

03 . 



Figure 3. The upper left panel shows the Local sample with SMBH mass plotted vs effective velocity dispersion (T02). The triangles and 
open circles represent SMBH masses determined via stellar dynamics and by other methods (masers, gas-disks, stellar proper motions) 
respectively. The grey symbols denote SMBHs whose mass estimates were deemed unreliable, or whose sphere of influence was not 
resolved (Marconi & Hunt 2003). Log-quadratic fits were performed on three samples. In the upper right panel we show a fit to the 24 
galaxies whose spheres of influence are resolved, and whose masses are deemed reliable. In the lower right panel we show a fit to the 20 
galaxies with SMBH masses determined via stellar dynamics. In the lower left panel we show a fit to 11 galaxies with masses determined 
from methods other than stellar dynamics. Projections of the contours represent the 1, 1.5, 2 and 2.5-sigma uncertainties on individual 
parameters. The 1-sigma error ellipse (thick line) for the full Local sample is shown superimposed on all three panels. The grey region 
shows the corresponding linear fit value for /3 from T02. In the upper left panel we show the three best fits. 


be wrong because they disagree with the a theoretical pre¬ 
conception that the Mbh — <J relation is log-linear. However 
concern has been expressed that that SMBH masses esti¬ 
mated using data that does not resolve the sphere of influ¬ 
ence are unreliable (e.g. Ferrarese & Ford 2004). Marconi 
& Hunt (2003) have ordered the SMBH sample in terms of 
the degree to which the SMBH sphere of influence is re¬ 
solved. There are 4 galaxies for which the sphere of influ¬ 
ence is smaller than the resolution (the half-width-at-half- 
maximum) of the data. In addition, Marconi & Hunt (2003) 
classed some of the mass estimates as unreliable for other 
reasons such as unknown disk inclination, resulting in the 
removal of a further 3 galaxies from the sample. In the up¬ 
per left panel of Figure 3 we show the remaining 24 galaxies 
(dark points). Those SMBH masses designated by Marconi 
& Hunt (2003) as being unreliable, and those with spheres 
of influence that are not resolved are plotted in grey. In the 
right panel of Figure 3 we show the log-quadratic fit to these 


24 galaxies. The result should be compared with the fit to 
the full sample, the 1-sigma error ellipse for which is also 
plotted (thick contour). We see that the removal of the 7 
unreliable galaxies has slightly increased the allowed range 
of /?2 • However the solution still prefers a non-linear term, 
the 1-sigma error ellipse lies near the /?2 = 0 line. 

In the lower right- and lower-left panels of Figure 3 we 
show the log-quadratic fits to the sample of galaxies with 
SMBH masses determined via stellar dynamics, and by other 
methods respectively. Each of these samples also shows evi¬ 
dence for a non-linear contribution to the A/bh — a relation. 
The sample containing SMBHs with resolved spheres of in¬ 
fluence, and the sample containing only SMBHs with masses 
determined via stellar dynamics have very similar solutions, 
indicating that SMBHs with masses determined via stellar 
dynamics in cases where the sphere of influence is unresolved 
do not skew the analysis and may be included in the sample. 

The question of the dependence of resolution on the 
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Figure 4. The upper left panel shows the Local sample with SMBH mass plotted vs central velocity dispersion (Ferrarese & Merritt 2000). 
The triangles and open circles represent SMBH masses determined via stellar dynamics and by other methods (masers, gas-disks, stellar 
proper motions) respectively. The grey symbols denote SMBHs whose mass estimates were deemed unreliable, or whose sphere of influence 
was not resolved (Marconi & Hunt 2003). Log-quadratic fits were performed on three samples. In the upper right panel we show a fit 
to the 24 galaxies whose spheres of influence are resolved, and whose masses are deemed reliable. In the lower right panel we show a fit 
to the 20 galaxies with SMBH masses determined via stellar dynamics. In the lower left panel we show a fit to 11 galaxies with masses 
determined from methods other than stellar dynamics. Projections of the contours represent the 1, 1.5, 2 and 2.5-sigma uncertainties on 
individual parameters. The 1-sigma error ellipse (thick line) for the full Local sample is shown superimposed on all three panels. In the 
upper left panel we show the three best fits. 


accuracy of the SMBH mass determined was investigated 
by Gebhardt et al. (2003). They re-evaluated the SMBH 
masses for galaxies where the SMBH mass had been previ¬ 
ously determined with ground-based imaging using higher 
resolution space-based data. Gebhardt et al. (2003) showed 
that i) the Magorrian et al. (1998) masses were indeed over¬ 
estimated, by up to a factor of 3, due to the assumption of 
2-integral models; and ii) comparison of 3-integral modeling 
of low-resolution ground based and high resolution space 
based data demonstrates that the SMBH mass determina¬ 
tions described in Gebhardt et al. (2003) contain no resolu¬ 
tion induced bias. In other words, provided that the model¬ 
ing is sufficiently general, the high resolution data increases 
the precision but not the accuracy of SMBH mass estimates 
(see figure 8 of Gebhardt et al. 2003). This implies that if 
the resolution were not sufficient, then the lower limit on 
BH mass would be < 0 (i.e. no definite detection). Thus the 
effect of spatial resolution is built into the estimate of the 


range of allowable mass. The masses of SMBHs in galaxies 
where the sphere of influence is not resolved should there¬ 
fore be just as reliable as those where it is resolved, in the 
sense that the determination is just as accurate but carries 
less precision (Richstone et al. 2004). 

At this point it should be noted that recent work from 
Valluri, Merritt & Emsellem (2004) has suggested that stel¬ 
lar dynamical estimates of SMBH masses are unreliable even 
when the SMBH sphere of influence is resolved. This unre¬ 
liability is due to a degeneracy between the mass-to-light 
ratio and SMBH mass, which was shown to become more 
prominent as the number of orbits used in the modeling is 
increased. On the other hand, Richstone et al. (2004) have 
countered this claim. They suggest that while it is true that 
the allowed range of SMBH mass increases with the number 
of orbits assumed, this range asymptotes to provide a true 
estimate of SMBH mass. Thus they conclude that provided 
the number of orbits is sufficient, SMBH mass estimation 
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from reconstruction of stellar orbits using kinematic data 
can lead to a robust SMBH mass estimate. Richstone et 
al. (2004) estimate the number of orbits that are required 
to achieve convergence in the mass estimate and conclude 
that published mass estimates have been made using a suf¬ 
ficiently large orbit library and are therefore reliable. 


Finally, the definition of the independent variable may 
be an important factor in the evaluation of any non-linear 
contribution to the Mbh — cr relation. A central velocity dis¬ 
persion has been advocated by Ferrarese & Merritt (2001), 
while an effective velocity dispersion was used by Gebhardt 
et al. (2000). The pros and cons of these different variables 
were summarised by Merritt & Ferrarese (2001). The central 
velocity dispersion is more easily measured than the effec¬ 
tive velocity dispersion which requires both surface photom¬ 
etry and spatially resolved spectroscopy over the effective 
radius of the galaxy. As a result, central velocity dispersions 
have been measured for many galaxies, while the effective 
velocity dispersions require significant effort to obtain the 
required data. On the other hand the effective velocity dis¬ 
persion better represents the velocity dispersion over the 
stellar spheroid, and may therefore be expected to offer a 
better representation of the depth of the spheroids gravita¬ 
tional potential well. Moreover T02 found that the SMBH 
itself could effect the value of the central velocity by up to 
30% in some cases. In order to evaluate the effect of changing 
the independent variable from an effective to a central veloc¬ 
ity dispersion, we repeat the above analysis on sub-samples 
of the Local sample of SMBHs (values of central velocity dis¬ 
persion are available in Merritt & Ferrarese 2001 & Ferrarese 
& Ford 2004). The results are shown in Figure 4, correspond 
to, and are qualitatively similar to to those in Figure 3, 
where the effective velocity dispersion was the independent 
variable. In particular, if the SMBHs whose spheres of in¬ 
fluence are not resolved are removed from the sample, then 
the sample still prefers a non-linear contribution. This non¬ 
linearity is at higher significance where the central rather 
than effective velocity dispersion is used. Similarly, samples 
of SMBHs with masses determined via stellar dynamics and 
by other methods each prefer a non-linear Mbh — o relation. 
This analysis shows the previously known result that use of 
the central velocity dispersion results in an estimate of the 
Mbh — cr relation that is steeper than where the effective 
velocity dispersion is used. 


In summary, samples where the SMBH masses were 
estimated by stellar dynamics, samples where the SMBH 
masses were determined via other methods, and samples 
where the SMBH masses are secure (Marconi & Hunt 2003), 
each have a fit that prefers the inclusion of a log-quadratic 
term. The sample of SMBH masses determined via other 
methods shows the strongest non-linear tendency, but with 
the largest uncertainty. This non-linear component of the 
Mbh —cr relation is present whether the independent variable 
is considered to be the central velocity dispersion (Ferrarese 
& Merritt 2000) or the effective velocity dispersion (T02). In 
light of these results, we consider only the full Local sample 
as described in T02, with the effective velocity dispersion as 
the independent variable, in the remainder of this paper. 


6 BAYESIAN PARAMETER ESTIMATION 


Our aim in this paper is to determine the parameters a, 
P and p 2 , as well as the intrinsic scatter (designated 8 ) of 
SMBHs around equation (4). Our general approach adds two 
features relative to the x 2 analysis described in the previous 
sections. First, our approach explicitly includes the intrinsic 
scatter (labeled 8 ) about the mean Mbh — cr relation as a 
third free parameter, and we consistently solve for all of 
a, p, p 2 and 8 . The second improvement is to include the 
asymmetric errors quoted for Mbh. 

In the absence of errors in Mbh and a we may define 
the contribution LHi of SMBH i to the likelihood LH at g^ 2t s 
for the set {a, p, P 2 ,8}. However both Mbh and a contain 
observational uncertainty. Each of the LHi must therefore 
be averaged over the uncertainty in both Mbh and a, hence 


LHi = 


d log (7 


dPi r 

dloga 


d log Mbh 


dPj 

d log Mbh 


LH{a, Mbh). 


( 6 ) 


The uncertainty in a is described by dPi/d\oga = 
N(d log at, 0.02) for all galaxies except the Milky-Way 
for which, following T02 we assume a larger uncertainty, 
dPi/d log cr = N(d. log <Ji, 0.08). Here the notation N(x\x,8x) 
refers to the value of a Gaussian distribution with mean x 
and variance Sx at x. The observed uncertainty in Mbh.i is 
assumed to follow a distribution 


dPj 

d log Mbh 


A(log Mbh| log Mbh.i, A log Mbh W i) 

if Mbh < Mbh,i 

A(log Mbh | log Mbh,i, A log M"PP) 

if Mbh > Mbh.i (7) 


Here we have defined A log M^? and A log Mbh)) as the un¬ 
certainty (in dex) above and below the observed value for 
SMBH i. The values of a, (3, @2 and 8 may then be esti¬ 
mated by maximising the product of likeli-lioods for the N g 
residuals 


LH, 


a,(3,@ 2 ,5 


= nSo LHi 


( 8 ) 


Note that we have not yet specified the definition of 
LH(a, Mbh). For a finite sample size, the solution for a, 
P, p 2 and 5 is sensitive to the definition of LH(a, Mbh). As 
a result, an inappropriate choice for LH(a, Mbh) can lead 
to a biased estimate of the parameters. In this paper we are 
analysing a sample where SMBH mass is unbiased at a fixed 
value of a. We therefore treat a as the independent variable 
and model the intrinsic scatter as a Gaussian with variance 
8 , hence 


LH(a, M bh ) = A(A Mbh ,i|0,5) 


where A Mbh) j = log M b h,< - 


a + P log 


+P 2 log 


200 km/s 

<n 

200 km/s 


(9) 


The joint a-posteriori probability distributions for the 
parameters a, p , P 2 and 8 may be found from 

d 4 P HP ■ HP ■ HP ■ HP ■ 

oc LH^'foj— -—-—-—, (10) 


dadpdP2d8 


dp dP2 


d8 
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Figure 5. The cumulative distribution of maximum likelihood solutions for the inferred values of /?i n fer and (32, infer in Monte-Carlo 
samples given different input values /?true and @ 2 ,true, which are shown by the vertical solid lines. The dots and dashed lines show the 
distribution means and medians respectively. Upper panels: The input relation had (3true = 4.2 and ,# 2 ,true = 1.6. The dashed lines show 
the 16th and 84th percentiles. Lower panels: The input relation had /3true = 4.02 and (32, true = 0. The dot-dashed lines show the values 
of (3 and (32 inferred from the data. 


where the P pr i or are prior probabilities, which we assume to 
be flat in this paper, i.e. 

HP • HP ■ HP • HP ■ 

'-*‘± prior '-*'± prior Li/.t prior prior /-i -i \ 

-^ oc ^r oc ^p oc ^^ ocl - (11) 

A-posteriori distributions for combinations of these param¬ 
eters may then be obtained by marginalising over the other 
dimensions, for example 


' dadpdfodS 


da J dt32 

1 */ —OO «/ —00 . 


dadfidfadS 


6.1 Bias in Parameter Estimation 


We generate samples of N s =31 SMBHs using the follow¬ 
ing procedure. For each of the N g = 31 galaxies we select 
a value of a drawn from the observed estimate N{cti, A ai). 
This value is used to select a SMBH mass from an input 
mean relation, offset randomly according to the input in¬ 
trinsic scatter <5 t rue- This mass is then further offset by a 
value drawn randomly from the quoted uncertainty in A/bh 
for the corresponding SMBH in the observed sample. We 
assumed input values of /3true = 4.2 and /? 2 ,true = 1.6. In the 
following section we show that these values lie near the best 
fit relation for the Local sample. The intrinsic scatter of the 
relations was assumed to be <5true = 0.3, defined in the y or 
ALbh direction. 

For each mock sample we find the most likely values for 
/3 and /?2 from equation (10). The cumulative distributions'!' 
of the best fit solutions are plotted in the upper panels of 
figure 5. The estimates are not skewed, but the parameters 


Before presenting the probability distributions for the pa¬ 
rameters /?, /?2 and S that describe the Local sample, we t Note that for a fair comparison, the most likely velocities from 

asses the bias in the fitting procedure by fitting parame- the sample are used in the fitting procedure, rather than the val- 

ters to Monte-Carlo realisations of mock SMBH samples. ues of a used to calculate the mock SMBH mass. 
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determined are slightly biased (by 8(3 = 0.1 and 5(32 = 0.2) 
as may be seen by the comparing the median (dashed line) 
and mean (denoted by the large dot) to the true value (solid 
vertical line) in each case. 


7 UNBIASED A-POSTERIORI PROBABILITY 

DISTRIBUTIONS FOR PARAMETERS IN 

THE LOG-QUADRATIC M BH cr RELATION 

Having assessed the bias in our parameter estimation, we are 
in a position to estimate the parameters describing the Lo¬ 
cal sample. In figure 6 we show the a-posteriori marginalised 
probability distributions for the parameters (3, 02 and 5 com¬ 
puted using equations (12-13) in combination with the like¬ 
lihood equation (9). The values of (3 and (32 corresponding 
to a fixed likelihood have been corrected a-posteriori for the 
biases A (3 = 0.1 and A/32 = 0.2. In the central row we 
show joint distributions for (32 and (3 (left) and for 8 and (3 
(right). The contours (dark lines) refer to 0.036, 0.14, 0.26 
and 0.64 of the peak height corresponding to the 4, 3, 2, 
and 1-sigma limits of a Gaussian distribution. In the bot¬ 
tom rows (dark lines) we show differential (solid lines) and 
cumulative (dashed lines) distributions for (3 (left) and 02 
(right). The vertical dotted lines show the variance. 

The contours show little correlation between parame¬ 
ters. Using the likelihood defined in equation (9) we find 
that the marginalised distributions imply (3 = 4.2 ± 0.37 
and 02 = 1.6 ± 1.3 (figure 6), and an intrinsic scatter 
of 8 = 0.28 ± 0.04. The normalisation was found to be 
a = 8.05 ± 0.06. The best fit (solid line) is plotted over 
the data in the upper left panel of figure 6, together with 
curves showing the level of intrinsic scatter (dotted lines). 
In the upper right panel of figure 6 we plot the residuals in 
log Mbh , together with horizontal dotted lines showing the 
value of the best fit intrinsic scatter. We have assumed an 
intrinsic scatter that is constant with a. Inspection of the 
residuals in figure 6 indicates that there is no systematic 
trend. 

The cumulative probability distribution for 02 shows 
that there is a positive contribution (at 90% confidence) 
from a log-quadratic term in the Mbh — cr relation that de¬ 
scribes the Local sample. This may be interpreted as indi¬ 
cating that at the 80% level the Mbh — a relation does not 
follow a single powerlaw between ~ 70km/s and ~ 380km/s. 
An alternative statistical question regarding the significance 
of this result concerns the frequency with which one might 
measure a log-quadratic term as large as the best fit of 
(3 2 = 1.6, assuming an intrinsically log-linear relation. We 
have performed fits to mock samples assuming an input in¬ 
trinsic slope of /3true = 4.02 and an input ,02, true = 0, which 
corresponds to the best fit log-linear relation (T02). The cu¬ 
mulative distributions for the best fit (3 and (32 are plotted 
in the lower panels of figure 5. We find that the estimate of 
(3 is unbiased in this case, and that the best fit value of (32 is 
greater than 1.6 around 10% of the time. This is consistent 
with the statement that the Local sample does not follow a 
single power-law at the ~ 80% level. 

The inclusion of a log-quadratic term in the fit leads 
to residuals that do not he systematically above the mean 
relation at high or low a. We can therefore use the be¬ 
haviour of these residuals to interpret the significance of 


the log-quadratic term. Our Monte-Carlo simulations (fig¬ 
ure 5) show that (after accounting for the bias) ~ 10% of 
samples with (3 = 4.2 and 02 = 1.6 produce best fits with 
(32 <0. As we saw in figure 1 the smallest three galaxies all 
have SMBHs that are more massive than the average im¬ 
plied by a log-linear fit. Each SMBH has ~ 1 chance in 2 of 
lying above the mean relation. This situation would there¬ 
fore arise by chance for ~ 1 sample in 10 (~ 1/2 3 ), which 
provides an intuitive understanding of the finding that there 
is a log-quadratic term in the Mbh — cr relation described by 
the Local sample that is positive at the 90% level. 


8 THE LOG-QUADRATIC ALbh — A/bulge 
RELATION 


In this section we apply our analysis to the sample of SMBHs 
and bulge masses summarised in Haering & Rix (2004). This 
sample closely resembles that described in T02. The SMBH- 
bulge mass relation is shown in the upper left panel of fig¬ 
ure 7, and can be described by a log-quadratic relation of 
the form 


log(Mbh) = Obulge + /3bulgelog(Mbulge/10 11 Af O ) 

+ (02, bulge [log(Mbh/lO n M 0 )] 2 . (14) 

We have repeated our x' 2 analysis on the Mbh — Mbuige sam¬ 
ple (not shown). We find that the slope estimated using a 
log-linear fit systematically increases by an amount that is 
greater than the statistical uncertainty when the smallest 
galaxies are removed from the sample. A x 2 minimisation 
of a log-quadratic relation provides a significantly better fit, 
with the power-law relation ruled out at the 1-sigma level. 
The parameters of the log-quadratic fit do not vary system¬ 
atically as galaxies are removed from the sample. 

Motivated by the success of the log-quadratic fit for the 
Mbh — Afbuige relation in a X' 2 minimisation, we perform a 
Bayesian analysis in analogy to the one described in the 
previous section. We define a likcli-hood function 


LHi = 


cL log Mbuige 


dPi 


d log Mbuige J _ 
dPi 


d log Mbh 


d log A4h 


LH (Mbuige, M bh ), (15) 


where LH(Mbu\ge, Mbh) is the likelihood of a set of A/bulge 
and Mbh given a relation described by abuige, 0 buige, 
/02, bulge and < 5 bui g e- Following Haering & Rix (2004) the 
uncertainty in Mbuige is described by dPi/d log A/bulge = 
N(d log ALbuige, 0.18) for all galaxies. At fixed bulge mass, 
we model the intrinsic scatter as a Gaussian with variance 
<5buige, hence we define a likeli-liood in analogy with equa¬ 
tion (9) 


LH (Mbuige, Mbh) = N(A Mbh ,i|0,<5buige) 
where A Mbh ,i = logM b h,i 

Albulge,* 


CTbulge T /3bulge log 


5 x 1O 1O M 0 


+02,bulge 


log 


Mi 


bulge, i 


5 x 10 10 AL o 


(16) 


We performed Monte-Carlo simulations as before and find 
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Figure 6. Upper Left: SMBH mass vs. central velocity dispersion from the compilation of T02. The solid line shows the best fit solution 
to equation (4) obtained using a likelihood defined in equation (9). The dotted lines show the level of intrinsic scatter around the best-fit 
relation. Upper Right: Residuals relative to the maximum likeli-hood solution. The dotted lines show the best fit value of intrinsic scatter. 
The lower panels show the a-posteriori probability distributions for the parameters (3, 02 and S. Center Left: Joint distribution for (3 and 
P 2 marginalised over a and S. Center Right: Joint distribution for S and (3 marginalised over a and (32- The contours refer to 0.036, 0.14, 
0.26 and 0.64 of the peak height corresponding to the 4, 3, 2, and 1-sigma limits of a Gaussian distribution. Bottom Left: Differential 
(solid line) and cumulative (dashed line) distributions for (3 marginalised over /? 2 , 6 and a. Bottom Right: Differential (solid line) and 
cumulative (dashed line) distributions for /?2 marginalised over (3 t S and a. The vertical dotted lines show the 16th, 50th and 84th 
percentiles. 
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Figure 7. Upper Left: SMBH mass vs. bulge mass from the compilation of Haering &; Rix (2004). The solid line shows the best fit 
solution to equation (14) obtained using a likelihood defined in equation (16). The dotted lines show the level of intrinsic scatter around 
the best-fit relation. Upper Right: Residuals relative to the maximum likeli-hood solution. The dotted lines show the best fit value of 
intrinsic scatter. The lower panels show the a-posteriori probability distributions for the parameters /?bulge> $2,bulge and ^bulge- Center 
Left: Joint distribution for 0hu\ge and /?2,bulge marginalised over c*buige an d ^bulge- Center Right: Joint distribution for ^bulge an d Amlge 
marginalised over Ubulge and 02, bulge- The contours refer to 0.036, 0.14, 0.26 and 0.64 of the peak height corresponding to the 4, 3, 
2, and 1-sigma limits of a Gaussian distribution. Bottom Left: Differential (solid line) and cumulative (dashed line) distributions for 
A>ulge marginalised over 02 , bulge? Amlge an d o;b u ige- Bottom Right: Differential (solid line) and cumulative (dashed line) distributions for 
02, bulge marginalised over 0, S and a. The vertical dotted lines show the 16th, 50th and 84th percentiles. 
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Figure 8. The values of logarithmic slope for the mean relations (solid lines) and for the relations using parameters on the edge of the 
1-sigma error ellipses that lead to the most extreme examples (dashed lines). The logarithmic slopes are shown as a function of a for the 
A/fah — a (left panel) and Mbh — A7 b ulge relations (right panel). The grey regions represent the 1-sigma ranges for ft (T02) and /3 b ulge 
(Haering & Rix 2004) using a log-linear relation. 


that the estimator described in equation (16) leads to 
a slightly biased estimate for the true values /3 b uige,true 
and fc, bulge,true- We find biases of A/Jbuige = 0.04 and 
A/?2, bulge = 0.035. 

In figure 7 we show the a-posteriori marginalised prob¬ 
ability distributions for the parameters /Ibuige, A,bulge and 
<5buige computed using the likelihood equation (16) and cor¬ 
rected for bias. In the central row we show joint distribu¬ 
tions for p 2 , bulge and Mbuige (left) and for <5 b uige and f3 hu i ge 
(right). The contours (dark lines) refer to 0.036, 0.14, 0.26 
and 0.64 of the peak height corresponding to the 4, 3, 2, 
and 1-sigma limits of a Gaussian distribution. In the bot¬ 
tom rows we show differential (solid lines) and cumula¬ 
tive (dashed lines) distributions for /3 b uige (left) and @2, bulge 
(right). The vertical dotted lines show the variance. We 
find that the marginalised distributions imply /3 b uige = 
1.15 ± 0.19, @2, bulge = 0.12 ± 0.14 and an intrinsic scatter of 
5 = 0.41 ± 0.07. The normalisation is a bu ige = 8.05 ± 0.1. 
The best fit relation has a positive log-quadratic term. How¬ 
ever this term is only required by the data at the 1-sigma 
level. Haering & Rix (2004) suggested an upper limit on the 
intrinsic scatter of ~ 0.3 dex, though the definition of this 
scatter was not specified. We have included intrinsic scatter 
self consistently in the analysis and find that the intrinsic 
scatter in SMBH mass at fixed bulge mass is <5 b uige ~ 0.4 
dex. This scatter is larger than the quoted uncertainties on 
either SMBH or bulge mass. The scatter in the M b h — M bu i g e 
relation is also ~50% larger than the scatter in the A/bh — cr 
relation. 

The best fit log-quadratic relation is plotted over the 
data (solid line) in the upper left panel of figure 7 together 
with curves showing the level of intrinsic scatter (dotted 
lines). In the upper right panel of figure 7 we plot the result¬ 
ing residuals, together with horizontal dotted lines showing 
the value of the best fit intrinsic scatter. Inspection of the 
residuals does not indicate any systematic trend of residuals 
as a function of M bu ige- 


9 DISCUSSION 

To quantify the extent of the departure of the Local sam¬ 
ple from a log-linear relation we have plotted the logarith¬ 
mic derivatives dlogM b h/dlog a and dlog ALbh/dlog ALbuige 
of the M b h — cr and Afbh — M bu ige relations as a function 
of cr and M bu ige in figure 8. The solid lines show the slopes 
of the best fit relations. The dashed lines show the slope of 
relations at the edges of the 1-sigma error ellipsoids in (3-/32 
(left panel) or /Jbuige — (32, bulge (right panel) space, while the 
grey stripes show the 1-sigma range for the slope of the log- 
linear fits from T02 and Haering & Rix (2004) respectively. 
The figure shows that in both cases the logarithmic slope of 
the log-quadratic relation varies substantially more over the 
range of cr or M bu i ge in the SMBH sample than the statisti¬ 
cal uncertainty in the slope (3 of a log-linear fit. This fact has 
contributed to statistical bias and inconsistent results when 
comparing the log-linear slopes of the M b h — cr relation be¬ 
tween different samples of SMBHs. In the remainder of this 
section we discuss various implications of a log-quadratic 
Mbh — u relation. 

9.1 The most massive SMBH 

The possibility of a log-quadratic M b h-cr relation leads to 
several interesting consequences. The first concerns the ve¬ 
locity dispersion of galaxies in which the most massive 
SMBHs reside. The mass of the largest SMBHs in the Sloan 
Digital Sky Survey (SDSS) survey volume was discussed in 
Wyithe & Loeb (2003). We repeat part of the discussion 
here. Consider a SMBH of mass M b h = 3 x 10 9 Mq. We 
can estimate the co-moving density of black-holes of this 
mass from the observed quasar luminosity function and an 
estimate of quasar lifetime. If we assume Eddington accre¬ 
tion, then at the peak of quasar activity, quasars powered 
by A/bh ^ 3 x 1O 9 AI0 SMBHs had a comoving density of 
4/ ~ 50Gpc -3 . Wyithe & Loeb (2003) found that the num¬ 
ber of quasars relative to the number of dark-matter ha¬ 
los at 2 ~ 3 implies a duty cycle for quasars of ~ 0.005. 
The comoving density of SMBHs with masses in excess of 
Mbh ~ 3 x 10 9 Mq at 2 ~ 3 was therefore ~ 10 4 Gpc~ 3 . As 
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Figure 9. Mass-functions for SMBHs. Two mass functions are 
shown. The grey stripe and dashed lines represent the ranges of 
uncertainty in the estimate of the mass-function assuming the log- 
linear and log-quadratic relations respectively. The ranges at fixed 
A/ hh in each case correspond to the 16th and 84th percentiles. The 
solid line shows the most likely estimate for the mass-function 
corresponding to the log-quadratic M^h — a relation. 


this corresponds to the peak of SMBH growth, the SMBH 
density should match the value observed today. This den¬ 
sity implies that the nearest SMBH of ~ 3 x 10 9 Mq should 
be at a distance dbh ~ 30Mpc which is comparable to the 
distance of M87, a galaxy known to possess a SMBH of this 
mass (Ford et al. 1994). What is the most massive SMBH 
that can be detected dynamically in the SDSS? The SDSS 
probes a volume of ~ lGpc 3 out to a distance ~ 30 times 
that of M87. At the peak of quasar activity, the density of 
the brightest quasars implies that there should be ~ 100 
SMBHs with masses greater than ~ 3 x 1 O 1 O M 0 per Gpc 3 , 
the nearest of which will be at a distance dbh ~ 130Mpc, 
or 5 times the distance to M87. The radius of gravitational 
influence of the SMBH scales as Mbh/cr 2 . We therefore find 
that for the nearest 3 x 10 9 Mq and 3 x 1 O 1 O M 0 SMBHs, 
the angular radius of influence should be similar. Thus the 
dynamical signature of the nearest ~ 3 x 1 O 1 O M 0 SMBHs 
on their stellar host should be detectable. 

In order to find these most massive SMBHs one needs to 
identify the galaxies in which they reside. We estimate the 
value of a implied for SMBHs of mass ~ 1 O 1 O M 0 using the 
mean Mbh — cr relation. If one adopts the mean log-linear re¬ 
lation (T02), then SMBHs with masses in excess of 1 O 1 O M 0 
should reside in galaxies with a <; 600km/s (see figure 11 ). 
However no such galaxies exist is the SDSS (e.g. Sheth et 
al. 2003), where largest values of galaxy velocity dispersion 
are found to be a ~ 400km s -1 . The mean log-quadratic re¬ 
lation reaches Mbh ~ 10 10 Mq at a more modest but still 
unobserved velocity dispersion, a ~ 500km/s. However the 
intrinsic scatter in the Mbh — & relation is <5 ~ 0.3 dex. We 
find that SMBHs of mass ~ 1 O 1 O M 0 differ by ^ 2(5 from the 
mean log-quadratic relation at a ~ 400km/s. The most mas¬ 
sive SMBHs with masses of ~ lO lo Af 0 inferred from quasars 
at z ~ 3 should therefore exist in the SDSS and would lie 
at ~ 25 above the extrapolated mean log-quadratic relation 
for galaxies with the largest measured velocity dispersions 
of a = 400km/s. This helps to reconcile the number of lumi¬ 


nous quasars observed at z ~ 3 with both the local Mbh — cr 
relation and the lack of super massive galaxies. 

9.2 The mass-function of SMBHs 

A positive log-quadratic term has a significant effect on 
the upper end of the SMBH mass-function. In figure 9 we 
show the SMBH mass-function computed by combining the 
Mbh — cr relation with the velocity dispersion function of 
Sheth et al. (2003). Two mass functions are shown. The 
grey stripe and dashed lines represent the ranges of uncer¬ 
tainty in the estimate of the mass-function obtained assum¬ 
ing log-linear and log-quadratic Mbh — cr relations respec¬ 
tively. The ranges at fixed Mbh in each case correspond to 
the 16th and 84th percentiles, and were computed taking 
into account Gaussian uncertainties in both the parameters 
of the Mbh — cr relation and the velocity-dispersion func¬ 
tion. The solid line shows the most likely estimate for the 
mass-function corresponding to the log-quadratic Mbh — cr 
relation. Figure 9 shows that inclusion of a log-quadratic 
term in the Mbh — cr relation results in an estimate for the 
number density of SMBHs with masses of ~ 1O 1O M0 that is 
several orders of magnitude larger than inferred from a log- 
linear relation. We have estimated the mass density (for a 
Hubbles constant of 71kms -1 Mpc -1 ) of SMBHs in the local 
universe. For a log-linear Mbh—cr relation we find a density of 
p = (2.30 ±0.45) x 1O 5 M0 Mpc -3 , while for a log-quadratic 
Mbh — cr relation we obtain p = (2.35±0.45) x 1O S M0 Mpc -3 . 
Hence the addition of a log quadratic term does not effect 
estimates of the total mass density which is dominated by 
SMBHs in the range 10' — 1O 9 M0 (Aller & Richstone 2002). 

9.3 Inclusion of SMBHs from reverberation 
mapping studies 

The sample of kinematically detected SMBHs will not grow 
by a large factor in the foreseeable future (Ferrarese 2003). 
Progress in understanding the statistical properties of the 
SMBH population may instead come via estimates of SMBH 
masses that are based on reverberation mapping studies 
(Gebhardt et al. 2000b; Ferrarese et al. 2001; Onken et 
al. 2004; Peterson et al. 2004; Nelson et al 2004). Onken 
et al. (2004) have presented a sample of 14 SMBHs with 
virial masses (A/bh.vir) determined via reverberation map¬ 
ping studies, in galaxies with measured velocity dispersions. 
Reverberation masses are determined up to a geometric fac¬ 
tor /, hence the SMBH mass is Mbh = f x Mbh.vir- Onken 
et al. (2004) have determined the average value (/) by com¬ 
paring the reverberation masses of these 14 SMBHs, to the 
masses of SMBHs in the Local sample using a log-linear 
Mbh —cr relation. They found (/) ~ 5.5 over a range of values 
for the slope /3 that include (5 = 4.58 (Ferrarese 2002) and 
/3 = 4.02 (T02). Assuming this scaling factor we can combine 
the 14 galaxies used to obtain (/) in Onken et al. (2004) with 
the Local sample of SMBHs to form a sample of 45 galaxies. 

The left panel of figure 10 shows the Mbh — cr relation 
for the Local sample (solid dots), and its comparison to the 
AGN SMBH masses and velocity dispersions (open circles) 
presented in Onken et al. (2004). We have repeated the \ 2 
minimisation using equation (4) for the combined sample 
of N g = 45 SMBHs. In the right-hand panel we show the 
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Figure 10. Left: The Local sample of SMBHs (solid dots) and the sample of SMBHs determined from reverberation mapping studies 
(with (/ v ) = 5.5; Onken et al. 2004; open circles). Right: The resulting contours of x 2 (minimised over a) of 0 and 02- Projections of 
the contours correspond to the 1, 1.5, 2 and 2.5-sigma limits on individual parameters. The point shows the most likely solution. 


1, 1.5, 2 and 2.5-sigma ellipsoids for 0 and 02 - The most 
likely solution has 0 = 4.05 and 02 = 1.8 with 1-sigma 
uncertainties of 0.3 and 1.0 respectively. The resulting best- 
fit Afbh — cr relation is over-plotted in the left panel (solid 
line). This solution is very close to the one obtained using 
only SMBHs in the Local sample (0 = 4.3 ± 0.35, 02 = 
1.9 ± 1.4). The intrinsic scatter required for a reduced x 2 of 
unity in this fit was eintrins = 0.25 dex, indicating that the 
addition of the AGN SMBHs does not increase the scatter of 
the Aibh — cr relation. The addition of AGN SMBHs results in 
a fit with a log-quadratic term which is non-zero at around 
the 1.5-sigma level. The sample of AGN SMBHs with masses 
determined from reverberation mapping therefore strongly 
supports a log-quadratic Mhh — cr relation. This is despite 
the fact that the normalization of the reverberation SMBH 
masses was determined via the assumption of a log-linear 
relation. 


9.4 The low mass Mbh — cr relation 

There are no SMBHs in the Local sample with masses below 
~ 2 x 10 6 ALq. Attempts to extend the Mbh — cr relation to 
lower masses through dynamical observations have centered 
on two cases, M33 (Merritt, Ferrarese & Joseph 2001; Geb- 
hardt et al. 2001), and NGC205 (Valluri, Ferrarese, Merritt 
& Joseph 2005). Only upper limits have been obtained in 
these systems. However we can compare the log-quadratic fit 
in the low-cr regime with SMBH masses obtained from single 
epoch virial estimates based on observations of AGN. The 
upper panel of figure 11 shows the extrapolation of the best 
fit log-quadratic relation for the Local sample, and its com¬ 
parison to the AGN SMBH masses and velocity dispersions 
presented in Barth, Greene & Ho (2004), as well as the point 
for POX52 (Barth, Ho, Rutledge & Sargent 2004). The AGN 
SMBH masses in this sample were estimated using the sin¬ 
gle epoch virial mass calibration of Onken et al. (2004). The 
log-quadratic fit derived using the Local sample appears to 
describe the SMBH masses in these AGN much better than 
the log-linear fit (see right-hand panel of Figure 11). The 
residuals are symmetrically placed around the best fit log- 
quadratic relation. On the other hand the dashed line in the 


right hand panel of figure 11 shows the additional residual of 
the log-quadratic relation relative to the log-linear relation. 
The residuals of the AGN SMBHs relative to the log-linear 
relation are systematically positive, indicating that the log- 
linear relation provides a poor description. 

The single epoch virial SMBH masses for the AGN 
have been estimated in a regime where the technique is not 
directly calibrated. One must therefore be cautious com¬ 
paring the AGN and local samples since the AGN SMBH 
mass estimates may contain some unknown systematic un¬ 
certainty relative to the dynamical estimates of the Local 
sample (Barth, Greene & Ho 2005 did not perform a fit to 
the combined sample for this reason). Moreover it is pos¬ 
sible that the AGN SMBH sample is biased towards large 
SMBH masses (Barth, Greene & Ho 2005). Never-the-less it 
is instructive to combine the Local and AGN SMBH samples 
and investigate the resulting log-quadratic A/bh — cr relation. 
We have repeated the \ 2 minimisation using equation (4) for 
the combined sample of N s = 47 SMBHs. In the lower-left 
panel of figure 11 we show the combined sample with the 
resulting best-fit Mbh — cr relation over-plotted (solid line). 
The intrinsic scatter required for a reduced \' 2 of unity in 
this fit was eintrins = 0.22 dex. In the lower-right panel of 
figure 11 we show the 1, 1.5, 2 and 2.5-sigma ellipsoids for 0 
and 02. The most likely solution has 0 = 4.25 and 02 = 1.3 
with 1-sigma uncertainties of 0.3 and 0.7 respectively. This 
solution is very close to the one obtained using only SMBHs 
in the Local sample (0 = 4.3 ± 0.35 , 02 = 1.9 ± 1.1 and 
eintrins = 0.25 dex), indicating that the AGN SMBHs fol¬ 
low the same log-quadratic ALbh — cr relation as the Local 
sample. However the addition of low-mass AGN SMBHs re¬ 
sults in a fit with a log-quadratic term which is non-zero at 
greater than the 2-sigma level. The sample of AGN SMBHs 
therefore support a log-quadratic ALbh — o relation. The log- 
quadratic relation describes the masses of SMBHs in galaxies 
with velocity dispersions ranging from 25-400kms _1 . If one 
extrapolates this to lower a then the log-quadratic relation 
offers the tantalising prediction that there is a minimum 
SMBH mass of ~ 10 s AL q, and that this minimum SMBH 
mass resides in bulges with a ~ 10km s -1 . Interestingly, 
this velocity dispersion corresponds to the minimum value 
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Figure 11. Upper Left: The best fit log-quadratic relation (determined from the Bayesian analysis of the Local sample) extrapolated to 
high and low cr (solid line), and the level of intrinsic scatter around the best-fit relation (dotted lines). Also shown (dashed line) is the 
log-linear relation from T02. The points from the Local sample are shown (solid dots). In addition, the open circles represent SMBHs 
with masses estimated using AGN identified in SDSS (Barth, Greene & Ho. 2004). The triangle represents POX52 (Barth, Ho, Rutledge 
& Sargent 2004). Upper Right: The corresponding residuals relative to the log-quadratic relation, with dotted lines showing the best 
fit intrinsic scatter. The dashed line shows the additional residual of the log-quadratic relation relative to the log-linear relation. Lower 
Left: The combined sample of Local and AGN SMBHs with the resulting best-fit log-quadratic relation determined from a x 2 analysis. 
Lower Right: The resulting contours of x 2 (minimised over of) of /3 and 02- Projections of the contours correspond to the 1, 1.5, 2 and 
2.5-sigma limits on individual parameters. The point shows the most likely solution. 


of the virial velocity for a dark-matter halo within which 
gas accreted from the IGM can cool via atomic transitions 
in hydrogen. 

9.5 The relation between 0 and /3buige 

Finally we would like to make the following point regard¬ 
ing the relationship between the power-law slope of the 
iV/bh — a and Mbh — A/b u i ge relations. An estimate of the 
bulge mass may be made using the virial mass, Mb u i ge ~ 
A/buige,vir oc cr 2 7? e , where R e is the scale radius of the 
bulge. The projection of the fundamental plane yields a re¬ 
lation between radius and velocity dispersion of R e oc cr 1,5 
(Bernardi et al. 2003). Consider the power-law slope at the 
mean velocity of the galaxy sample (a ~ 200km/s). We have 
Mhh oc (Mbuige)^ 1 " 11 ® 6 , which implies Mbh oc u 3 ' 5 ^>ui g G or 
0 = 3.5/3buige- This is in excellent agreement with the best 
fit parameters of 0 ~ 4.2 and /Jbuige = 1.15 derived in this 


paper. Furthermore, a value of 0 = 3.5 for the Mbh — cr 
relation is ruled out at ~ 80% significance. If it is to be 
consistent with the A/bh — cr relation, the A/bh — Mbuige rela¬ 
tion should therefore be steeper than linear for galaxies with 
cr ~ 200km/s. Given the relation between cr and the virial 
mass, we would also expect the A/bh —A/bulge relation to have 
a positive quadratic component, given a positive quadratic 
component in the Mbh — cr relation. 


10 CONCLUSION 

Using the Local sample of SMBHs we have demonstrated 
that the parameter describing a log-linear fit to the A/bh — cr 
relation is sensitive to the inclusion of individual galaxies 
at a level larger than the statistical uncertainty. This indi¬ 
cates that a log-linear A/bh — cr relation does not provide a 
good description of the Local sample of SMBHs. We expand 
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a general relation between log(Afbh) and log(cr/ 200 kms _1 ) 
to second order, and fit the data using the resulting log- 
quadratic relation instead. We find that a log-quadratic re¬ 
lation provides a substantially better fit to the Local sample 
of SMBH masses and velocity dispersions. Moreover unlike 
the log-linear relation, the parameters describing the log- 
quadratic relation are not systematically dependent on the 
inclusion of individual galaxies in the sample. 

After allowing for a second-order term in the Mbh — cr 
relation we find an unbiased estimate for the slope of the Lo¬ 
cal sample at a = 200km/s to be (3 = 4.2 ± 0.37. This value 
is slightly (2/3-sigma) larger than previous estimates for this 
sample. However the logarithmic slope dlog A/bh/dlog a of 
the best fit log-quadratic relation varies substantially, from 
2.7-5.1 over the velocity range of the Local sample. The coef¬ 
ficient of the second order term (32 = 1.6± 1.3 is greater than 
zero at the 90% level. This indicates that with 80% confi¬ 
dence the Local sample describes an A/bh — a relation that 
does not follow a single powerlaw between ~ 70km/s and 
~ 380km/s. We have tested the sensitivity of this conclu¬ 
sion to different sub-samples of SMBH masses determined 
via different techniques, as well as SMBHs with and without 
resolved spheres of influence. The log-quadratic fit implies 
a non-linear contribution to the Mbh — & relation in a sub¬ 
sample of galaxies that contain only SMBHs whose spheres 
of influence were resolved, in a sub-sample containing only 
SMBHs with masses determined through stellar dynamics, 
and in a sub-sample containing only SMBHs with masses 
determined only by non stellar dynamical methods. More¬ 
over the non-linearity is present in each sub-sample whether 
the central or effective velocity dispersion is used as the in¬ 
dependent variable. 

SMBH masses in active galaxies estimated via reverber¬ 
ation mapping offer an avenue to increase the SMBH sample 
for study of the Mbh — o’ relation. We find that the combina¬ 
tion of the 14 galaxies with reverberation SMBH masses and 
measured velocity dispersions (Onken et al. 2004) with the 
Local SMBH sample leads to a log-quadratic relation with 
the same best fit as the Local sample alone. 

The A/bh — o relation can be extended to lower masses 
through the inclusion of single epoch virial estimates of 
SMBH masses based on observations of AGN (Barth, Greene 
& Ho 2004). In a \ 2 analysis the best-fit log-quadratic 
Mbh — o’ relation is unchanged by the addition of a sam¬ 
ple of 16 low mass SMBHs. However the uncertainty in /3 2 is 
reduced. We find that in the absence of a systematic error in 
the normalisation of SMBH masses between the Local and 
AGN samples, the A/bh — o relation described by the com¬ 
bined sample deviates from a power-law at greater than the 
2-sigma level. The best-fit log-quadratic relation predicts a 
minimum mass for SMBHs in galaxies of ~ 10 s Mg, which 
should reside in bulges with a ~ 10 km s -1 . 

A log-quadratic A/bh — cr relation has important impli¬ 
cations for SMBH demography. In particular, estimates of 
the local SMBH mass-function that utilise the log-quadratic 
Mbh — c relation describe densities of SMBHs with M b h 
lO 9 A /0 that are orders of magnitude larger than expected 
for a log-linear Mbh — cr relation. In addition the departure 
from a power-law should provide important clues regarding 
the astrophysics responsible for the Mbh — cr relation. For 
example one recent model describing the effects of radiative 
feedback on SMBH growth predicts departure from a power- 


law relation, including steepening of the relation at large cr 
(Saznov et al. 2005). 

We have also applied our analysis to the relation be¬ 
tween SMBH and bulge mass using the sample described in 
Haering & Rix (2004). We find evidence for a log-quadratic 
term in the Mbh — A/bui ge relation (at the 1 -sigma level), with 
A>uige = 1.15 ± 0.19 and /3 2 ,bulge = 0.12 ± 0.14. We find an 
intrinsic scatter of dbuige = 0.41 ±0.07 dex. Since we find the 
intrinsic scatter in the Mbh — cr relation to be J = 0.28±0.04 
dex, there is ~ 50% more scatter in the SMBH mass at fixed 
bulge mass than at fixed velocity dispersion. 

The sample of kinematically detected SMBHs will not 
grow by a large factor in the foreseeable future (Fer- 
rarese 2003). Progress in understanding the statistical prop¬ 
erties of the SMBH population will instead come via es¬ 
timates of SMBH masses that are based on reverberation 
mapping studies. The increased sample of SMBHs will offer 
the possibility of more clearly defining the local Mbh — cr re¬ 
lation (Gebhardt et al. 2000b; Ferrarese et al. 2001; Onken 
et al. 2004; Peterson et al. 2004; Nelson et al 2004), and of 
extending its study to high redshift (Shields et al. 2003). 
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APPENDIX A: COMPARISON OF 
PARAMETERS FOR DIFFERENT SMBH 
SAMPLES 

Following the discovery of the A/bh — cr relation (Gebhardt 
et al. 2000, Ferrarese & Merritt 2000), there has been de¬ 
bate in the literature among the different groups regarding 
inconsistencies in the value of the measured slope. These 
differences have been attributed to several causes including 
use of biased statistics (Merritt & Ferrarese 2001) and sys¬ 
tematic differences in the velocity dispersions used for the 
host galaxies (T02). We have shown that the slope (3 of a 
log-linear relation is substantially changed by the presence 
or absence of one or a couple of the smallest galaxies in the 
sample. Given the different samples of velocity dispersions 
that have been employed by the different groups, the bias 
introduced by the assumption of a log-linear relation may 
therefore be more or less significant for different samples. 

Merritt & Ferrarese (2001) presented an analysis of the 
log-linear relation for three different samples. The first two 
samples were drawn from Ferrarese & Merritt (2000, sam¬ 
ple 1), and from the additional galaxies presented by Geb¬ 
hardt et al. (2000, sample 2). The accuracy of some of the 
additional SMBH masses included by Gebhardt et al. (2000) 
has been called into question on the basis that the SMBH 
sphere of influence was not resolved (e.g. Merritt & Fer¬ 
rarese 2001). In addition, Merritt & Ferrarese (2001) dis¬ 
cuss the combined sample from these two sub-samples. For 
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Figure Al. Log-quadratic fits to the samples of Ferrarese & Merritt (2000, sample 1), the additional galaxies added by Gebhardt et 
al. (2000, sample 2), as well as a combined sample (upper right panel) as described in Merritt & Ferrarese (2001). Projections of the 
contours represent the 1, 1.5, 2 and 2.5-sigma uncertainties on individual parameters. The 1-sigma error ellipse for the combined sample 
is repeated in the lower panels. In this figure the correlation is between SMBH mass and the central velocity dispersion (Ferrarese & 
Merritt 2000). The grey region shows the corresponding linear fit value for 6 from Merritt & Ferrarese (2001). The upper left panel 
shows the three best fits along with the data. Sample-1 is represented by open circles, and sample-2 by triangles. The large cross shows 
the point j3 = 4.9, 62 = 3.75 which lies inside the 1-sigma error ellipses of all three samples in both the cases where central and effective 
velocity dispersion are used as the independent variable. 


sample 2 Merritt & Ferrarese (2001) computed aperture cor¬ 
rected central velocity dispersions, so that all galaxies in the 
combined sample could be considered when estimating cor¬ 
relations between SMBH mass and both a central (Ferrarese 
& Merritt 2000) and an effective velocity dispersion (Geb¬ 
hardt et al. 2000). Merritt & Ferrarese (2001) present fits for 
the slope j3 of a log-linear relation to these samples using re¬ 
gression with bi-variate errors and intrinsic scatter (their 
label BRS). 

We first consider the correlation of SMBH mass with 
central velocity dispersion (Ferrarese & Merritt 2000). Mer¬ 
ritt & Ferrarese (2001) find /3 — 4.81 ± 0.55 and (3 = 
3.75 ± 0.59 for samples 1 and 2 respectively. For the com¬ 
bined sample they get f3 = 4.43±0.39. Samples 1 and 2 differ 
at the 2-sigma level, while the combined sample leads to a 
slope that lies between the two. To mimic this analysis for a 
log-quadratic relation we minimise the following y 2 statistic 


2 _ y- {yi-a- fixi - /3 2 x 2 ) 2 
* Z -4 ltrins +eT + (/3 2 + 2 [ > 

where yi and Xi are the logarithm of SMBH mass in solar 
masses and the logarithm of velocity dispersion in units of 
200km/s respectively. The variables e x i and e y i are the un¬ 
certainties in dex for these parameters. The variable eintrins 
is the intrinsic scatter, adjusted to yield a reduced y of 
unity for the best-fit solution. In figure Al we show the 1, 
1.5, 2 and 2.5-sigma error ellipsoids for (3 and (3 2 in sample-1, 
sample-2 and the combined sample. The 1-sigma ellipsoid for 
the combined sample is shown in the other panels for com¬ 
parison, and the best fit relations are plotted over the data 
in the upper left panel. The grey regions represent the 1- 
sigma range for (3 determined by Merritt & Ferrarese (2001). 
The ellipsoids show that all three samples suggest a posi¬ 
tive quadratic term is present in the Mbh — a relation. The 
values of f3 determined by Merritt & Ferrarese (2001) are 
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Figure A2. Log-quadratic fits to the samples of Ferrarese & Merritt (2000, sample 1), the additional galaxies added by Gebhardt et 
al. (2000, sample 2), as well as a combined sample (upper right panel) as described in Merritt & Ferrarese (2001). Projections of the 
contours represent the 1. 1.5, 2 and 2.5-sigma uncertainties on individual parameters. The 1-sigma error ellipse for the combined sample 
is repeated in the lower panels. In this figure the correlation is between SMBH mass and the effective velocity dispersion (Gebhardt et 
al. 2000). The grey region shows the corresponding linear fit value for (3 from Merritt & Ferrarese (2001). The upper left panel shows the 
three best fits along with the data. Sample-1 is represented by open circles, and sample-2 by triangles. The large cross shows the point 
(3 = 4.9, /?2 = 3.75 which lies inside the 1-sigma error ellipses of all three samples in both the cases where central and effective velocity 
dispersion are used as the independent variable. 


conditional probabilities for (3 given (32 = 0, and so should 
correspond to regions where the error ellipsoids cross the 
/?2 = 0 line. Figures Al and A2 show this to be the case. 
The inclusion of a quadratic term in the relation leads to 
similar estimates for parameters using each sample, with 
(3^5 and (32 ~ 4 lying inside the 1-sigma ellipsoids. 

We also investigate the correlation between SMBH mass 
and effective velocity dispersion (Gebhardt et al. 2000). 
The values for effective velocity dispersion a in Gebhardt 
et al. (2000) did not include an estimate of the uncer¬ 
tainty. In the absence of quoted uncertainties, Merritt & Fer¬ 
rarese (2001) assumed different values for the uncertainty in 
effective velocity dispersion and showed that increasing the 
uncertainty lead to steeper best fit slopes, and to consistency 
between fits for samples 1 and 2. T02 has since advocated 
5% for the fractional uncertainty in effective velocity disper¬ 
sion. We assume 5% errors in a and adopt a = lOOkm/s for 
the Milky-Way. In this case Merritt & Ferrarese (2000) found 
(3 = 3.6±0.5 and (3 = 4.11±0.35 for samples 1 and 2 respec¬ 


tively. For the combined sample, they found f3 = 4.05 ± 0.3. 
In figure A2 we show the results for the log-quadratic fit 
in this case. The results are consistent for a log-quadratic 
relation, with (3 ~ 4.6 and (32 ~ 3 lying inside the 1-sigma 
error ellipses of all three samples. Consistency between the 
best fit log-quadratic relations for samples 1 and 2 therefore 
does not require uncertainties in a to be larger than 5%, the 
value advocated by the Nuker team (Gebhardt et al. 2003). 

Finally we compare the estimates of (3 and (3-2 for re¬ 
lations between SMBH mass, and the central or effective 
velocity dispersions using the combined sample. The use of 
an effective velocity dispersion (with 5% uncertainty) leads 
to smaller values of f3 and (32 relative to a fit using central ve¬ 
locity dispersions. Merritt & Ferrarese (2001) found that this 
difference could be removed by increasing the uncertainty in 
a to greater than 10%. However a solution with (3 ~ 4.9 and 
P2 ~ 3.5 lies inside the 1-sigma error ellipsoids in both cases. 
Overall we find that a solution with (3 ~ 4.75 and (32 ~ 3.7 
(shown by the large cross in figures Al and A2) lies on or 
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inside the 1-sigma error ellipsoids of all three samples using 
either the central or the effective velocity dispersions as the 
independent variable. Therefore if a log-quadratic Mbh — & 
relation is used we find that there is no significant disparity 
between the different samples. 


APPENDIX B: AN UPDATED SAMPLE OF 
SMBH MASS AND CENTRAL VELOCITY 
DISPERSION 

In a recent review, Ferrarese & Ford (2004) have compiled 
an updated sample of local SMBHs, listed in their table-II. 
Of the 38 mass estimates, they class 8 as unreliable. Of the 
remaining 30 SMBHs, the mass observations resolved the 
sphere of influence in only 25 cases. Ferrarese & Ford (2004) 
used this sample of 25 reliable mass estimates where the 
sphere of influence is resolved to generate a log-linear fit to 
the Afbh — cr relation. They find a log-linear slope of (3 = 
4.86 ± 0.43. We have performed log-quadratic fits to this 
updated sample, and show the results in Figure Bl. In the 
upper left panel we show the sample of 25 SMBHs with 
resolved spheres of influence, and the corresponding best- 
fit relation. In the upper right-hand panel we show the 1, 
1.5, 2 and 2.5-sigma error ellipsoids of f3 and 02 for this 
sample. This sample prefers a log-quadratic relation, with 
the log-linear relation lying near the A\' 2 = 1 contour. The 
grey regions in Figure Bl represent the 1-sigma range for 
/3 determined by Ferrarese & Ford (2004). We see that the 
contours of \ 2 intercept the = 0 line over a range of slope 
described by the log-linear solution. 

Of the 5 galaxies with reliable mass-estimates, but un¬ 
resolved spheres of influence in Table II of Ferrarese & 
Ford (2004), 4 have had their stellar kinematical data mod¬ 
eled using both high-resolution spaced based data, and lower 
resolution ground based data (Gebhardt et al. 2003). Geb- 
hardt et al. (2003) found that while higher resolution imag¬ 
ing led to a more precise determination, there was no sys¬ 
tematic trend of the best fit values. We therefore repeat the 
above analysis with the inclusion of the 5 galaxies and show 
the results in the lower row of Figure Bl. The additional 
5 galaxies only change the slope f3 of a log-linear fit by a 
small amount (significantly less than the error bar). More¬ 
over, these 5 SMBH masses are evenly spread about the 
best-fit relation, indicating that there is no systematic bias 
due to their spheres of influence not being resolved. However 
the inclusion of these galaxies significantly increases the sig¬ 
nificance of a log-quadratic term in the Mbh — cr relation. 
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Figure Bl. Upper left : The sample of 25 SMBHs with reliable masses and resolved spheres of influence described in Ferrarese & 
Ford (2004). Upper right: The resulting contours of (minimised over a) of (3 and /? 2 - Projections of the contours correspond to the 1, 
1.5, 2 and 2.5-sigma limits on individual parameters. The point shows the most likely solution. Lower left and right: As for the upper 
row, but including the additional 5 SMBH masses in table II of Ferrarese & Ford (2004) classed as reliable but with unresolved spheres 
of influence (shown as dark points). The grey regions show the allowed range of (3 = 4.86 =t 0.43 in a log-linear fit to the 25 SMBHs with 
a resolved sphere of influence (Ferrarese & Ford 2004). 
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